Damped collective motion of many body systems: 
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Abstract 

We address the problem of collective motion across a barrier like en- 
countered in fission. A formula for the quantal decay rate is derived which 
bases on a recently developed variational approach for functional integrals. 
This formula can be applied to low temperatures that have not been acces- 



fN| I sible within the former PSPA type approach. To account for damping of 

Q!^ ■ collective motion one particle Green functions are dressed with appropriate 

^ . self-energies. 
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1 Introduction 



o , _ 

^ , For a particle moving in a metastable potential J.S. Langer [Tj has shown that 

the decay rate 

i? = MM^ (1) 



><: -T 

can be calc 
free energy 



' can be calculated from the partition function Z via the imaginary part of the 



J^=-T\nZ. (2) 

(As usual temperature is measured in MeV such that the Boltzmann constant is 
put equal to = 1.) The parameter 1/k (where k < 0) sets the time scale for 
the solution {q)t = {q)t=o exp(— Kt) which describes the motion of the average 
trajectory away from the barrier top. Using path integrals Langer's formula 
has extensively been exploited in the framework of the Caldeira-Leggett model 
(CLM) |2I for dissipative quantum systems. Within this model the decay rate has 
been studied in great detail from high to very low temperatures. The region of 
the "crossover temperature" Tq, below which quantum tunneling becomes more 
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important than pure thermal activation has been treated in OHIIH]; for reviews 
on this subject see e.g. [HllZl- 

In the present paper we want to adapt similar methods to systems like finite 
nuclei, where the dynamics of the particles is governed by a varying mean field. 
Different to the CLM this field, and hence the "heat bath", changes with the 
variable which parameterizes the path along which the free energy exhibits the 
barrier. A first attempt along these lines was reported in . There, the partition 
function was calculated by accounting for correlations around the mean trajectory 
which were of RPA-type. In different context such an approach to Z has been 
phrased "Perturbed Static Path Approximation" (PSPA), for details see [HI HHl 
ITT] . As this method is based on a quadratic expansion of the effective action about 
stationary points it suffers the same breakdown as an analogous expansion within 
the CLM; the rate diverges at the Tq. In the present paper we like to overcome this 
problem by employing a variational approach to the quantum partition function, 
suggested recently in [T21 HHI • 

At finite temperature nuclear collective motion is damped. In [8^ this feature 
was simulated by applying simple energy smearing to the results obtained within 
the independent particle model (IPM). In the present paper we like to apply a 
more microscopic picture. The latter bases on the replacement of the one-body 
Green function of the IPM by dressed ones, accounting in this way for incoherent 
scattering processes among the particles. In this way connection is established 
to transport theory for damped collective motion (see e.g. [EllTn]). Indeed, by 
using linear response functions we are able to introduce microscopically transport 
coefficients like the inertia M, friction 7 or the stiffness C, which then may serve 
as input for the calculation of the decay rate of metastable systems. 

The paper is organized as follows: In Sect. El we briefly review several approxi- 
mations to the partition function Z of interacting many body systems. Thereafter 
in Sect.Elwe take into account residual interactions, apply linear response theory 
to treat damped harmonic motion and establish the connection to microscopic 
transport theories. In Sect. El we modify the methods of Sect. El such that dis- 
sipation can be treated. After the derivation of formulas for the decay rate of 
metastable interacting many body systems in Sect. El we want to estimate the 
accuracy of our results. In the ideal case this should be done with an exactly 
solvable model for the many body system. However, to the best of our knowledge 
for damped collective motion, the case we are interested in, such a model does 
not exist. For this reason so far we are only able to test in Sect. IHlour approach 
at the example of a particle moving in a one-dimensional potential under the 
influence of damping, leaving a more involved analysis for many body systems 
for future work. Finally we like to close by discussing our flndings and drawing 
some conclusions. 
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2 Approximations to the partition function 



The collective degree of freedom shall be introduced through a mean field ap- 
proximation to a separable interaction of the type 16^ kF'^/2. It is to be added 
to the unperturbed Hamiltonian H, which at first shall be assumed to be of one 
body nature. The total 7i thus is given by 

n = H+^FF. (3) 

For k < the interaction is attractive and leads to iso-scalar modes in the nuclear 
case. The operator F is meant to represent an effective generator of collective 
motion along a fission path. For the evaluation of the partition function Z we 
want to exploit functional integrals in imaginary time t = . . .hjS with (3=1 /T. 
The two body interaction may then be treated by the Hubbard-Stratonovich 
transformation (HST) to the mean field Hamiltonian 

nns^[q{r)]=H + q{T)F. (4) 

Here, the g(r) represents the collective variable. Following ^\ it will be treated 
by an expansion around an average value go defined through hf5qQ = J^^ dTq{T). 
The deviation will be expanded into the Fourier series ^(t)— go = J2rjto 1r exp(iz/rr), 
with the Matsubara frequencies = (27r//i/3) r and q-r = ql- As described else- 
where (see e.g. [Hll CHI 113) the partition function can be expressed by the 
following (ordinary) integral over the go: 



Z{P) = y ^ / "^^0 exp(-/3^^P^(/5, go)) C(/5, go) (5) 

In this expression jF^^^ is the go-dependent effective free energy in the quasi-static 
picture, left over here in the classical limit, where all r-dependent deviations from 
go can be neglected. For convenience we introduce the effective free energy 

^^*^(/3, go) = ^^P^(/3, go) - ^ lnC(/?, go) , (6) 

which includes quantum effects through the factor C,. This correction still is a 
functional integral in gj.-space and specified through an Euclidean action se(/3, go). 
Various existing approximations consist in the way this action is treated. As 
already indicated through the notation employed in (0) and ©, within the Static 
Path Approximation (SPA) the C^^^ is put equal to unity, which implies to discard 
all fluctuations around go. Corrections to this approximation (see e.g. [12]) are 
obtained by expanding systematically the action 

Prst QrQsQt + ^ CTrstuqrqsqtqu] + 0{q^) (7) 
' ' \r-,S7^0 r,s,t=/=0 r,s,t,u^O / 
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around the static path go to some power in the g^, exhibited here to fourth order. 

Within the SPA+RPA |9', the Perturbed Static Path Approximation (PSPA) 
[TU] or the Correlated Static Path Approximation (CSPA) 11^ the expansion ((7j) 
is truncated after the second order 



^e'^^ = ^E^-^^^- (8) 



This approximation takes into account quantum fluctuations around go on the 
level of local RPA modes. It turns out (see e.g. |HI) that the second order 
coefficients of the PSPA are diagonal A^s = Aj./2 Sr-s, with the being given 
by 

Here, the response function x{^) has been introduced, which is defined through 
the relation 

6{F)^ = -xH5qiLu) (10) 

where Sq{uj) means a general time-dependent deviation from go. Mind that the 
averages have to be calculated from the Hamiltonian at fixed go. In this 
way the x{^) depends on /? and go. Since this x{^) oiily contains excitations of 
nucleonic nature, 

hhJki = ek-ei, (11) 

which are to be distinguished from collective ones (to which we will come to be- 
low), we may call it the nucleonic or intrinsic response function. In the following 
the single particle energies shall be measured with respect to the chemical poten- 
tial fj,, which means to write hjJk{qo) = €k{qo) — l^'i.Q.o)- The local RPA frequencies 
zUfj,{P,qo) appearing in ^ must be calculated from the secular equation 

l + kx{zu,) = 0. (12) 

Note that the important factors in the product are those which deviate from 
one. This happens whenever the w"^ is sizably different from the corresponding 
uli- This is so in particular for the typical collective modes, for which for stable 
isoscalar modes, for instance, one has zu"^ <^ uli- 

For proper convergence of the path integral for (^pspa needs the condition 
(see e.g. [HI CHI) 

Ai(/?,go)>0. (13) 

It ensures the PSPA to be well defined such that the quantum corrections can be 
written in the form 
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The restriction ()13|) can easily be translated into a condition on that temperature 
at which an unphysical divergence of the quantum fluctuations around go occurs: 
T > To. Here 

^|tX7inst(go)| 

2n 

can be calculated from the unstable local RPA mode for which ^f^^^^ < 0. The 
latter is present whenever J^^^^{qo) develops a barrier. 



To = Max<j„ — (15) 



Obviously, the unphysical divergence at To is a deficiency of the PSPA. It can 
be cured by accounting for inharmonicities in the expansion ((7j) of the action. To 
be able to lower the breakdown temperature to To/2 it is sufficient to concentrate 
on the three coefficients pi,i,-2, P-1,-1,2 and cti, 1,-1,-1 and neglect the remaining 
ones JH|- This improved approximation was called extended Perturbed Static 
Path Approximation (ePSPA). For T > To/2 the quantum correction factor may 
then be expressed by the form 

C'^'^P, qo) = exp [x\P, go)] erfc [x{P, go)] H j;^^ ■ (16) 

The quantity 



^1 4^ ^ " 6ai,i,-i,_i - ^' ^ (17) 

is given in terms of the third and fourth order inharmonic terms of the Euclidean 
action ((Zj). The coefficients p and a can be expressed through one body Green's 
functions 



h Z — UOk 

corresponding to the static part of the Hamiltonian Details of the calculation 
can be found in ^Sl^l- With Fermi occupation numbers denoted by n{tk) the 
general coefficient arstu consists of a sum of terms like 



\k\ 

IT 



i,k,m,o 



+ n(efc) gf\ujk + iVs) gTi'^k + i«^r+s) g^i'^k - m) 
The prst or the higher order coefficients are of similar structure. 



2 APPROXIMATIONS TO THE PARTITION FUNCTION 



6 



Recently for the many body system studied here, a variational approach to 
the quantum corrections C of (0) has been developed in [T21 CSl- Based on the 
Euclidean action of the PSPA jHI) a reference action 



^^Y.^r\qr\'. (20) 



\k 



r>0 



is introduced, where the RPA frequencies zuf^ appearing in Q are replaced by 
variational parameters 

The are adjusted such that the expectation value (se — sf^f^ with respect to 
the Sq and consequently the effective free energy © is minimized. Accounting in 
the action Q for inharmonic terms up to fourth order and using the abbreviation 

within this approximation the correction factor of the variational approach can 
be written as [El Cnj 

In C(^) (/?, qo) = -Y. 1^^^ - E - y E r r • ^^^^ 

T->0 r>0 r,s>0 ^ 

Here, the fourth order coefficient arstu again is calculated according to (fT^ . For 
convergence of the double sum in the second line of (j23|l the behavior of (Trstu 
at large r, s, . . . is crucial. By inspection of (jlHI), using the unperturbed Green's 
functions (fTK|) and z/^ ~ r, one can convince one-self that arstu falls off rapidly 
enough to assure convergence. 

This variational approach IT^ to the partition function of an interacting 
many body system improves the accuracy of the PSPA and the ePSPA. Moreover, 
even in cases where the static free energy !F^^^ develops a barrier (as a function 
of go), this novel approach still is applicable below To or To/2. Evidently, this 
variational approach is the easier to handle the fewer variational parameters 
one needs to introduce. For the present formulation the latter are unequivocally 
related to the RPA frequencies of the secular equation (fT^ . In principle, the 
latter has as many solutions as there are p-h excitations. However, as already 
discussed below dH)), not all of them are equally important for collective motion. 
In the ideal case the latter is dominated by only a few or eventually by just one 
prominent mode. Then it suffices to concentrate on these or this latter one(s) 
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3 Accounting for residual interactions 



So far the operator H appearing in (jH} and is identified with that of the 
independent particle model (IPM), with the corresponding excitations given by 
(fTTj) and the one body Green's functions by (fTHj). The response function defined 
by (fTUj) can be written in the form x{^) = x'(^) + ^x"{^)- Calculated from the 
IPM its imaginary part x" consists of a sum of (5-functions located at discrete 
energies. However, the IPM neglects the residual two body interaction between 
the nucleons, which describes the incoherent scattering of particles and holes and 
consequently couples the Iplh excitations of to npnh excitations and true 
compound states. It is this residual interaction which in the end is responsible 
for genuine relaxation processes and, hence, for damping of collective motion. 
To account for such effects we will take a pragmatic point of view. Rather than 
to embark on the difficult (if at all feasible) task to account for the residual 
interaction in full glory we will simply take over the factorization of the response 
function into one body Green's functions gk{z), which using the corresponding 
spectral densities gk{z) leads to 

— n{nn + ^) [gk{z + n)Qi{Q) + Qk{n)gi{n - z)] . 

(24) 

For a more detailed description we refer to ^3] where also references to earlier 
work can be found. Different to the IPM we dress the Green's functions by 
self-energies = — irfc/2 and write 

^'^^'^ = lz-uj,-Mz)/h- ^^^^ 
instead of (fTHj). These self-energies are calculated from 

^k{z) = / 26 

i_oo 2-^ z-n 

with the following phenomenological ansatz for the widths 

As for the two parameters introduced this way, the values Pq ~ 33 MeV and 
c ~ 20 MeV have been used in the past (see e.g. [Hj and [IS]). It can be said 
that at zero temperature the widths calculated this way are in good agreement 
with empirical data for single particle excitations ^H] • Moreover, at T = forms 
of this type have been used in [201 in analyses of experimental results within 
optical model type approaches, for more details see Sect. 4.2.3. of [T3]. 
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We should like to note that the factorization assumption leading to (|2^ is 
applied to describe the intrinsic dynamics but not that of collective modes. One 
may therefore argue that coherent effects are small, in particular at larger ther- 
mal excitations. After all the underlying SPA is meant to represent the high 
temperature limit. 

We may proceed now to establish connection to microscopic transport theories 
like the one described in [T3]. To this end it us useful to introduce the collective 
response function defined by 

5(F)^ = -XcollM/extH, (28) 

where fext stands for an external field coupled to via fext{T)F. As is well 
known (see e.g. [E] or ^^1) it can be obtained from the nucleonic response 
function of ()10|) through 

Its poles are given by the solutions of the secular equation (fT^ . Very similar to 
the situation discussed in connection with Q, not all those poles cjj^ contribute 
equally strongly. In the ideal case one of them is shifted down to frequencies which 
are much smaller than the corresponding nucleonic excitations (|tJ7;xP ^ ooli for 
stable iso-scalar modes). Moreover, it may acquire a considerable part of the 
overall strength. 

In order to take into account dissipative effects we perform two modifications. 
As mentioned previously, the nucleonic response function x{^) is calculated from 
the dressed Green's functions fl^ instead of (fTHjl . Using the T of ^T7\i the imag- 
inary part {x^)"{uj) becomes a continuous function of u. As we are interested 
in a description of slow iso-scalar collective motion, as given for instance for nu- 
clear fission, we will concentrate on the lowest lying (pair of) poles of the (locally 
defined) collective response function (jSHI)- This truncated form is fitted by the 
response function of a damped harmonic oscillator of inertia M, stiffness C and 
damping 7 

-1 _ 1 -1 

~ Mu^ + i-fu-C ~ ]C\ {uj/wy + 2ir]uo/w - sgnC ' 

Then, instead of several ro^ only one collective mode is left over. The frequency 
w and the effective damping coefficient t] are defined as follows: 

^ = ITT and 71 = — — — = , (61) 

M ' 2Mw 2^/M\C\ 

The parameter rj indicates underdamped {j] < 1) or overdamped motion (77 > 1). 
After these modifications, which may be summarized as 

xZf{^) — ^ Xcoii(^) ^ Xosc(^) , (32) 
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the secular equation (fT^ for the local collective modes reduces to the equation 

(-Y + 2ir](-) -sgnC = 0. (33) 

Its solutions are given by the two frequencies 

uj± = tu ^±A/sgnC — r/^ — . (34) 

In contrast to the RPA frequencies zu^ of the IPM, the uj± have a finite imaginary 
part. This implies a damped local collective motion due to the influence of the 
residual interaction. 

At this place it is important to stress that the transport coefficients C = 
C{T, go), ^ = viT, go) and w = w{T, go) or M = M(T, go) and 7 = 7(T, go) have 
been derived from a microscopic quantum theory starting from the separable two 
body interaction (jS)). Like the response function itself, all of them depend on 
the temperature T = and the collective coordinate go- With respect to 
nuclear collective motion this dependence may be considered a great advantage 
over the Caldeira-Leggett model (CLM) P|, which is often used to describe open 
quantum systems (see e.g. |7j; for a criticism from nuclear physics point of view, 
see [inilll])- There, the transport coefficients are not calculated microscopically. 
Rather they are introduced as parameters which are independent of T and go. 
This latter feature goes along with the fact that in the CLM the set of intrinsic 
degrees of freedom, which acts as the "heat bath" for the collective ones, is 
modeled by a set of fixed harmonic oscillators. The latter do not vary with the 
collective degree of freedom, and so does the spectral density. Such an assumption 
is not applicable to nuclear fission. There the heat bath for collective motion is 
given by the nucleonic degrees of freedom which move in a mean field which 
depends on the nuclear shape and thus varies with the collective variable. 



Dissipative systems: Modification of the ap- 
proximations to the partition function 



In this section we study how the inclusion of dissipative effects modifies the 
approximations put together in Sect. El in formal sense. First we look at the 
PSPA. After the reduction (jH^ to a single damped mode an evaluation of 

1 

\r = l + k x{}Vv) = ; 7^ (35) 

(see (0) and (1221)) leads to 

A. - XP = + ^g"^,^' . (36) 
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Here, the (inverse) time scale 
and the quantity 



— = T7 (37) 



C — k \h\ 
w\P,qo) = -j-r- = sgnCzu\P,qo) + , , ' ' . > sgnCzu\P,qo) (38) 
M M[p,qQ) 

have been introduced. A calculation of the nucleonic response function 

= Mu^ + {C - k) 

from ()30p by inversion of the formula (j^^ shows that w'^ plays the role of the 
nucleonic frequencies. In dissipative PSPA the form of the quantum corrections 
to the classical partition function defined in (jS)) does not change as compared to 
the IPM (HD). Only the are replaced by the A^.'^^ of The convergence 

condition (fT^ now turns to A^^''(/3, go) > and explicitly reads 

f2TiV 1 27r _ . 

\n(3j Tkin h(3 

In the case of stable modes with C > this condition is always fulfilled. For 
unstable modes with C < 0, in analogy to (fTH|l the local crossover temperature 
To (go) is given by the smallest root of the l.h.s. of the inequality pn|l . Using the 
effective damping rj of pip this To (go) becomes 

To(go) = ( V^TTTM - Viq.)) ■ (41) 

Its maximal value To = MaXg(,To(go) defines the global crossover temperature. 
Like in the CLM for dissipative quantum systems (see e.g. [SUHIIZ!) we find that 
damping r/ > diminishes this To. 

Concerning the anharmonic terms of the Euclidean action there is an important 
difference between our approach and the CLM. In the latter the conservative 
and the dissipative forces are decoupled from each other. As one implication, 
for instance, the inharmonic terms of the potential V'^"'\q) with n > 2 are not 
modified by dissipative effects. In the case of interacting many body systems 
the situation is more complicated: The replacement of the unperturbed Green's 
functions (fTHjl by the dressed ones (pSj) not only leads to a finite damping strength 
7] and a change in the frequency of the local harmonic motion uj±. Due to the 
structure (fT^ of the expansion coefficients of the Euclidean action ((7j) also the 
inharmonic terms are influenced by this replacement, such that p^st ~^ Prst 

(T) 

<^rstu <^rstu- Taking into account dissipative effects for the quantum corrections 
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of the ePSPA along the hnes of one can easily rederive (fTBj) . but with — 
X^^ and 

^ - - /it • («) 

where 

i^^^^ = 6ai2-,-. - . (43) 

Using the variational approach in order to account for dissipative features the 
form 

..o_M^A(r)|^^|2 (44) 



\k\ 



r>0 



with 



is used as reference action. The trial action ()44p contains only one variational 
parameter VL^ and, analogously to the case sketched in Sect. emerges from the 
dissipative version of s^^^ (see (jHl) with X^P of (j36|l instead of A^) by the re- 
placement sgnC Ci7^ VL^. Using and repeating the derivation of 13] for the 
dissipative case the quantum correction factor to the classical partition function 
in the variational approach turns out to be given by (j^Hj) with the replacements 

and (Jrstu cri^tu everywhere: 

In C^^^ go) = - E InAp -J^U^-^^Yl ^^l-s (47) 

r>0 r>0 ^ r,s>0 ^« 



Please note that this form is very similar to that of the extension of the 
Feynman-Kleinert variational approach (FKV) |2I] to open quantum systems 
(see e.g. |7]). There the quantum partition function of a particle of mass M in a 
one-dimensional potential V{q) and an oscillator heat bath environment has been 
studied within the CLM. Using the (inverse) time scale ^7\i and the fluctuation 
width [Zj 

2 . 1 

= E ^2 + ^JrU^,^.) + n2^f3, q,) ^^^^ 
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to fourth order in a the quantum corrections read [7]: 

In C^^^ (A go) = 2^ / f''\^n2m ^ (49) 

[V"{qo)/M - go)] a\P, go) 

-^r(^)(go) (a2(Ago))^ + 0(a6) 

RecaUing the definitions of and Ilr^'' in P5|l and (j46|) respectively, the terms 
of (jTfjl and correspond to each other in the order of their appearence. In 
(jTfjl the analog of the squared width pH|l is a sum consisting of terms of the form 
1/A)- ' (see (gni)), which converge to 1 in the limit r ^ oo. As mentioned already 
before, different to the fourth order term in ()49|) for convergence of the double 
sum in the last line of (jlTj) the behavior of cr^^l^.^ for large r and s is crucial. 
This difference can finally be traced back to the different path integral measures 
in both cases. 



5 Decay of metastable states 

In jS] a formula for the decay rate of damped interacting many body systems 
has been derived within the PSPA. Like that whole formalism it is applicable 
for P < Po only. For P ^ Pq the quantum corrections to the classical rate show 
an unphysical divergence. From the fact that the variational approach to the 
partition function can be applied also for /3 > /5o we expect that a rate formula 
can be derived, which behaves well in the crossover region and beyond. 

Based on hanger's findings PP for the decay rate i? of a metastable system in 
the literature (see e.g. jH] or [7j) for high temperatures the formula 

i?(/?</?o) = -||-Im^(/3), (50) 
n po 

is in wide use. With the identifications 1/Po = Tq (form (PT|1 evaluated at the 
barrier go = g^) and \k\ = —iu^ (from (jMjl with w = zub and sgnC = —1) it 
becomes identical to (0). For the low temperature region I. Affleck has shown 
[22] that 

R{P>Po) = -'^lmJ^{P) (51) 

gives the same result as a Boltzmann average over the energy dependent decay 
rate R = —2/hlmE. Of course, the formulas (jSU)) and coincide at /? = Pq. 
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In the following we apply (jSUj) and ()51|) to our modeling of interacting many 
body systems by the Hamiltonian To this end the integral (0) for the partition 
function and its relation to the free energy of the total system (j2)) is used. This 
integral is dominated by the free energy J^SP^(go) of the SPA. In the sequel the 
latter will be assumed to have just one minimum ai q = qa and one barrier at g = 
qb- (Whenever suitable the indices a and b will be used for quantities which have 
to be evaluated at the minimum and the barrier, respectively.) Moreover, we want 
to concentrate on examples where the barrier is sufficiently pronounced. First of 
all this implies that the (temperature-dependent) height B^^^ = J^^^^ — J-'^^^ 
is sufficiently large, as compared to temperature, j3B^^^{j3) 3> 1. The minimum 
and the barrier are assumed to be well separated in go- In addition we assume 
that locally they may be approximated by oscillators with stiffnesses C^^^ = 
d"^ J-"^^^ / dq^ , with a positive C^^^ and a negative C^^^. Any further structure of 
j^sPA jg neglected. Moreover, the variation of the quantum corrections lnC(/3, go) 
with go should be sufficiently weak such that the integral © can be evaluated by 
a steepest descent approximation. The partition function 

Z{(3)^Z,{P) + iZb{P) (52) 

then develops an imaginary part which is associated to the stationary point at 
the saddle and given by 

= exp[-f3:Fr\(3)] UP) . (53) 



The real part 



comes from the stationary point at the minimum. To obtain a special 
treatment is necessary for the integration around the barrier top. Following 
Langer P], there the contour of integration must be distorted into the upper 
complex half plane. Because of the assumptions about the barrier height we have 
lmZ{i3) <C IieZ{(3), implying that the imaginary part of the free energy can be 
approximated by 

lmT{(3) ^ -1/13 ■ ImZ(/?) /ReZ(/?) . (55) 

Using (jKn|) or (jKTjl the rate can be separated into a product of two factors, one 
representing the classical limit and the other the quantum corrections: 

i?(/3) = /2dass(/5) ■ /qm(/3) (56) 



5 DECAY OF METASTABLE STATES 



14 



Inserting the crossover temperature Tq = l//3o of ()4H) and the imaginary part 
of the free energy into (fSUj). for high temperatures the classical part reads: 



(57) 

This result has aheady been found in 0. There it also had been rewritten in a 
more intuitive way by making use of (j3T| for the stiffnesses: 




In this way in (j37j) besides the Arrhenius factor exp[—PB^^^] the attempt fre- 
quency Wa appears, with which the system tries to overcome the barrier. Notice 
please that in this rate formula the factor in brackets is equivalent to Kramers' fa- 
mous correction factor fl^. It indicates that the classical contribution to the rate 
is reduced by damping. Different to the rate formulas derived in the CLM in 
the ratio of the inertias at the barrier and the minimum appears jH] • It should be 
mentioned, however, that this additional factor is not seen in hanger's original 
work although his derivation started from a Fokker-Planck equation where the 
inertia is allowed to depend on the collective variable. To clarify this point one 
may eventually have to modify the way the collective degree of freedom is intro- 
duced in the common SPA and PSPA formulations of the HST on which also our 
approach bases heavily. 

For low temperatures application of leads to 



RtE^siP >(^o) = ^^ exp[-/3i?sP^(/3)] . (59) 

Note that this region has not been accessible with the PSPA in 0. The formulas 
(I57|) and ()59|1 describe the contribution of thermal activation to the decay rate. 

Finally, we turn to the quantum correction factor /qm(/5). Both for (jKTj) as well 
as for it has the same form, 

Mm - ^ (60) 

and increases the rate. In [H] it has already been pointed out for the PSPA 
that using (|Tl|l with (j36|l the factor (j60|l has better convergence properties than 
the ad-hoc generalization of the CLM to coordinate-dependent damping. In the 
following for the evaluation of ( and /qm we will limit ourselves to the variational 
approach. 
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6 An estimate for the accuracy of the rate for- 
mula 

In Sect. El the framework for the calculation of the decay rate of metastable self- 
bound Fermi systems such as heavy atomic nuclei has been developed. To test 
the accuracy of the formula (jSU)) with ^7\i or and (jSHl) it would be desirable 
to have at one's disposal a model for metastable many body systems, for which an 
exact result can be given even in the case of finite damping. Unfortunately, to the 
best of our knowledge, such a model does not exist. One might think to take the 
Lipkin-Meshkov-Glick model (LMGM) j;24j in the version of |lSj. Indeed, in this 
model the free energy J-'^^^{qo) develops a barrier below a critical temperature. 
This barrier separates two minima in a symmetric, bound system. Therefore it is 
not possible to introduce a decay rate in proper sense. Also based on the LMGM 
in [22] an exactly solvable model for many body systems has been developed, 
which allows to study the fission process. Unfortunately, both models do not 
allow for exact solutions at finite damping. 

Because of this problem we will follow a different strategy and look at the 
much simpler one-dimensional case of a particle of mass M moving in the cubic 
potential 



with minimum at qa = 0, barrier at qb = 2Q/3 and with a barrier height B = 
Vb—Va = 2 Mzu^ Q'^/27, although our approach is suitable for many body systems. 
For this situation the quantum decay rate has been calculated exactly with help 
of the CLM for all temperatures T = . . . oo and various damping strengths |21 El 
mE]. In addition, in the FKV there exists an analog of the variational approach 
used in the present paper, namely for a particle moving in a one-dimensional 
potential. 

Unfortunately, in its original form [21] the FKV is not directly applicable to 
the potential (jfilj) . The reason is that in the expansion ()49|1 only derivatives of 
V{qQ) of even order contribute. Consequently, like for the harmonic oscillator, for 
a cubic potential the FKV reduces to a Gaussian approximation and inharmonic 
terms are not seen. To cure this problem, in [2^] besides the trial frequency Q a 
second variational parameter has been introduced, which controls the minimum 
of the local oscillator. Here, we follow a simpler strategy and consider the stable 
fourth order potential 



which by adjusting n > 1 and C can be fixed such that it reproduces the barrier 




(61) 




(62) 
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region of ()61|) to arbitrarily high accuracy. By apphcation of the formulas ()5Up 
and ()51|). for which the free energy has to be evaluated for the particle of mass 
M in the potential (j62|) . one gets (see e.g. jHlIZI) 



J\Wl'\/M / ; \ / W" 

i?ciass(/3 < /3o) = 2n (v^^+^-^) y exp{-l3B) (63) 



RciUP>(3o) = ^y^exp(-/35) (64) 

instead of (|K7jl - (fHnj) . Using (jlHI) for the calculation of the correction factor (fHlljl 
one can get an estimate on the accuracy of the formulas derived in Sect. El 




Figure 1: The metastable potentials in one dimension, for a barrier height of 
B = ?)hw. The coordinate q is given in units of Q and energies are scaled to 
hw. The effective potential iy^^^(/5, q) is shown for [3 = 2(3o and two effective 
damping strengths rj = and rj = 1. 



In Fig. ^we show the barrier region of the potentials V{q) and W{q,n = 3) 
with barrier height B = 3huj. The global minimum of the potential W{q, n = 3) 
is found at q^ ~ 2.4(5 i'^ot shown in the figure) and is of value —W{q+) ~ 
A2hw ^ B. Because of the large depth of this minimum a thermally activated 
flux originating from this minimum is strongly suppressed. Also shown in Fig. ^ 
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is the "effective classical potential" j^I] 



W^^^^(/?,go) = W{qo) - - lnC^^^(/?,go) (65) 

for two effective damping strengths and a given temperature. It is the analog of 
(jni) evaluated in the variational approach and contains quantum corrections C^^^? 
which are taken from (jl^ with V replaced by W. As compared to W{qo) the 
effective classical potential W^^^{(3, go) has a smaller barrier and slightly smaller 
stiffnesses at qa and qb- 




Figure 2: A comparison of different approaches to the decay rate for a particle 
in a one-dimensional cubic potential with barrier, for zero [r] = 0) and critical 
[f] = 1) damping. The arrows point to the /3o = I/Tq where the PSPA breaks 
down. For more details see text. 



In Fig. 121 we show the /3-dependence of the rate formulas in logarithmic rep- 
resentation, for two effective damping strengths. The results of two different 
calculations are presented. The one denoted by "FKVl" (fully drawn line) is 
evaluated by using formulas (jU^ and with the quantum corrections /qm 

taken from ()60p and PU)) . For the second calculation ("FKV2", fully drawn line 
with symbols) the saddle point approximation is applied directly to the effective 
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classical potential (j^ . The quantal rate is calculated from the formulas 



FKV 



1 / (vrJKV)" 



where B^^^ = W^^^ — W^^^ is the (temperature-dependent) effective barrier 
height within the FKV, which is smaller than B, as can be seen from Fig. ^ 
Please note that FKV2 is a less drastic approximation to the go-iiitegral for Z 
than FKVl. In addition we show in Fig. |21 the classical rate as a dotted line, 
for which f^^^ = 1. The analog of the PSPA rate for a particle moving in 
a one-dimensional potential takes into account quantum effects via a Gaussian 
approximation. In order to keep the figure transparent these curves are not shown 
here. The inverse crossover temperature (3q = I/Tq, however, where the PSPA 
formalism breaks down due to a divergence of f^^^, is marked by arrows. For low 
temperatures P > Pq also the result of H. Grabert et al. jH EI is plotted (dashed 
line). Taking from Tab. I and II of ^ the relevant quantities, these results were 
obtained from an exact numerical treatment of the dynamical "bounce solution" 

The classical decay rate shows the purely exponential behavior known from Ar- 
rhenius' law. The rates derived from FKVl and FKV2 coincide with the classical 
result at high temperatures (small /3) but show enhancement at larger /3. In the 
crossover region /5 ^ /5o these approximations smoothly connect the classical rate 
with the result of Grabert et al. for P > Pq. For very low temperatures (/5 ^ Pq), 
however, FKVl and FKV2 deliver bad results, too. Instead of converging to a 
finite limit, both versions qualitatively behave more like the classical rate and fall 
off exponentially. This can be understood as follows: Also the FKV is basically 
a static approximation that incorporates inharmonic terms by smearing out the 
potential locally, see Fig. ^ This works well as long as the quantum fluctua- 
tions are not too large, as given for moderate to high temperatures. At very low 
temperatures, however, there is no way around accounting for the full nonlin- 
earity of the dynamical bounce solution qB{T), which minimizes the Euclidean 
action for P > Pq and therefore gives the leading contribution to the decay rate. 
Due to this limitation of the variational approach, it would simply be asking too 
much to expect good agreement with the fully quantum rate even at extremely 
small temperatures. As the FKV2 contains more information about the quantum 
effects than the FKVl it underestimates the rate less dramatically at low tem- 
peratures. Especially in the region /3 > /?o FKV2 still delivers acceptable results 
where FKVl is no longer reliable. 
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